Dynamical suppression of telegraph and 1/ / noise due to quantum bistable fluctuators 
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We study dynamical decoupling of a qubit from non gaussian quantum noise due to discrete 
sources, as bistable fluctuators and 1/f noise. We obtain analytic and numerical results for generic 
operating point. For very large pulse frequency, where dynamic decoupling compensates deco- 
herence, we found universal behavior. At intermediate frequencies noise can be compensated or 
enhanced, depending on the nature of the fluctuators and on the operating point. Our technique 
can be applied to a larger class of non-gaussian environments. 
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Controlling the dynamics of a complex quantum sys- 
tem is at the heart of quantum information [lj . However 
in any real device the computational variables entangle 
with the environment leading to decoherence Q. Bang- 
Bang (BB) control techniques have been proposed as a 
way to achieve an effective decoupling from the environ- 
ment 0, 0]. They may be operated by a sequence of 
strong external pulses separated by a time At 0. For 
At — > full decoupling 0,|j] of the unwanted interactions 
is achieved. The physics in this limit is a manifestation 
of the quantum Zeno effect ■ 

In practice At is finite especially when full-power 
pulses are used. This imperfect decoupling is still well 
described by the Zeno limit if At <C 7 _1 , the typical 
time scale of the environment . If 7 is large one may 
argue that BB chops noise and frequencies Li < 1/ At < 7 
are averaged out. This optimistic scenario could foresee 
applications to solid state coherent devices, where low- 
frequency noise la] is the major problem for quantum 
state processing pi 0] • Investigation of this point is 
one of the topics of this Letter, where we study environ- 
ments of dissipative quantum bistable fluctuators [jj . 

Recently decoupling from classical Random Telegraph 
Noise (RTN) was studied in the Zeno limit, At <C 
7 _1 Gaussian noise with 1/f spectrum has also 

been studied ^lj and decoupling for decreasing At was 
found. On the other hand in echo protocols, details of 
the structure and of the dynamics of a solid state discrete 
environment [f3 | may become important if the condition 
At <C 7 _1 is not met. 

We consider a qubit (TLq = — | a z — — a x ) coupled to 
an impurity. The Hamiltonian is 

H = H Q - X -g z E + He + V(t) (1) 

The environment Hamiltonian He = Hd + Ht + Hb 
describes an impurity level occupied by a localized elec- 
tron, TLd — e c Sd, tunneling with amplitudes (/Ht = 
SfcTfccj[d + h.c.) to a fermionic band, described by 
Hb = Ylk £ k c k c k- The charge in the impurity is cou- 
pled to the qubit, E = vSd. Control is operated as in 



Ref. 0|, the external field V(t) being a sequence of ir- 
pulses about x. This model may describe charge noise 
due impurities close to a solid-state qubit 0, 13 0] • 
The characteristic scale of the impurity is the switching 
rate 7 = 27rA/"(e c )|T| 2 (Af is the density of states of the 
fermionic band, |T fe | 2 » \T\ 2 ). 

This environment is non gaussian a key feature to 
explain recent observations in Josephson qubits (split- 
ting of spectroscopic peaks, beats in the coherent oscilla- 
tions [lfij ) due to individual impurities close to the device. 
The observed 1/f noise is due to a set of such impuri- 
ties |9j- We find that decoupling of this environment is 
sensitive to details of its dynamics. If pulses are not very 
frequent it shows a rich variety of behaviors, suggesting 
that BB may also be used for spectroscopy. 

We operate with instantaneous pulses which do not 
modify the environment, the corresponding evolution op- 
erator being Sp ~ ia x (g) He- The evolution operator of 
the Hamiltonian is [SpS] 2N , where S — exp(—iH.At) 
with V(f) = is the evolution between pulses. Echo cor- 
responds to N = 1. The reduced density matrix (RDM) 
of the qubit is obtained by tracing out the environment 

p(t) = Tr E {[S P S] 2N W(Q) 1&S P ] 2N } = £tW)\ (2) 

where W(t) is the full density matrix and £t[-] is the quan- 
tum map [l| associated to the reduced dynamics starting 
from a factorized state, W(0) — p(0) ® w E \ll\ . 

We may try to app roximate Eq.(J2J) by a Bloch-Redfield 
master equation ■ In this framework the environment 
remains in equilibrium and the map for the RDM in the 
first At has the Lindblad form £a*[p(0)] ~ exp(£A£)p(0). 
The factorized structure of W(t) is preserved if we apply 
pulses, so subsequent At can be treated in the same way. 
After t = 2NAt we get 

p(t) « [V ■ e^] 2N p(0) (3) 

where V is the superoperator of the pulses. This approx- 
imation, which is correct for a weakly coupled and fast 
environment, yields that BB does not affect the decay 
of the coherence. Of course BB decoupling is effective 
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only in situations where memory effects are paramount, 
and the trace in Eq.Q must be taken at the end of the 
protocol. In this cases we should go beyond the approxi- 
mation Eq.©. The possibility we explore is to treat part 
of the environment on the same footing of the system • 
We denote with p{t) the RDM of the qubit plus localized 
level. The system is now described by Hq — \o~ z E + Hd 
and we use the same steps leading to Eq.©. The ma 
p(t+At) = e CAt p(t) is evaluated by a master equation 
Ht being the interaction and Hb the bath. The RDM of 
the qubit p Q = Tr<z[p(i)] is obtained by tracing out the 
localized level at the end of the protocol. 

We express p(t) in the basis \Q^) = \8 n ±)\n), where 
\n) (n = 0, 1) are eigenstates of d)d and \8 n ±) are the 
two eigenstates of Hq — (v/2)na z , their energy splitting 
being fi„ = y/(e + nv) 2 + A 2 . We denote |o) = |6»+), 
1 6) = \0„), \c) = \9+) and \d) = |0f ). In Refs. it was 
found that the impurity remains in an unpolarized state, 
TiQ{p(t)} = J2 n=0 iPn(t)\n)(n\, if initially this was the 
case. This simplifies the dynamics of pif. the only non 
vanishing entries are the four populations and the coher- 
ences Pab(t) and p c d{t) (with the conjugates). Thus we 
should diagonalize a 8 x 8 sub-matrix of C. Using the rep- 
resentation of V, this is enough to find the approximate 
map Eq.(PJ) for a BB protocol, at all times. 

If A = 0, the calculation can be carried out analyti- 
cally. In absence of pulses [H, <J Z ] = 0, the populations 
of the qubit do not relax while its coherences are given 
by (0 O + |Tr d [p(t)]|&o-> - Pab(t) + p cd (t). This holds 
true also for an even number of pulses. This symmetry 
further simplifies C leading to independent evolutions of 
populations (subscript p) and coherences ((ft) 
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where C p /<f, are 4x4 matrices, whereas is a 2 x 2 
matrix acting on the vector p^ = (p a b, Pcd)- The pulse V 
is also diagonal in the p — i> indexes. In the 0-subspace 
it is given by / ® o~ x , which allows to obtain the map for 
coherences (p^, p^*) in an echo procedure 
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The "diagonal form" implies that the game reduces to the 
two-component map p${t) — [e r * At e r, * At ] P4>(0). This 
can be cast in a convenient form if the map expiT^t) 
found in Refs. Q is represented in SU(2) 



^(t) = [D/|a| 2 



(5) 



Here a = [(1 — iw) 2 — 2iSp eq g — g 2 } 1 ^ 2 determines the 
rates of the multi-exponential reduced dynamics of the 
qubit, the parameter g — (f2i — ^o)/7 quantifies non- 
gaussianity 0, 5p eq is the equilibrium population differ- 
ence of the fluctuator and w is related to the energy shifts 




FIG. 1: For a fixed t = IO7"" 1 we plot T N (i)/g 2 for BB proce- 
dures with TV echo pair of pulses. The parameter is a = v/y. 
N — corresponds to Free Induction Decay (FID) [al^. A 
gaussian environment with the same power spectrum would 
give, for arbitrary g, the curve here labeled with g = 0.1, 
since Tjv(t) oc p 2 Inset: Fjv(t) for g = 1.1 for different At 
(lines with dots, 7Ai = 5, 2, 0.2) are compared with the FID 
To(t) (thick line) and with results obtained by a stochastic 
Schrodinger equation (crosses) simulated with a very efficient 
piecewise deterministic algorithm |2j, Il9j| . Simulations are not 
accurate at relatively long times and in general they require 
a large statistics for the process (we used 10 6 realizations). 



produced by the band. Finally T>(At) = (T> x ,V y ,'D z ) 
and the quantities X>j(At) are easily found from the re- 
sults in Refs. (eg. V (At) = |a| 2 |cosh^| 2 + (1 + 
g 2 + w 2 )|sinh2^*| 2 ) and determine D(At) = [V - 
PI 2 ] 1 / 2 , x(At) = arctanh(|P|/I? ) and ovg = a ■ V/\V\. 

Eq. (J5J allows to discuss the dynamic decoupling of a 
quantum bistable fluctuator. We expect a rich physics 
since this environment has distinctive features depend- 
ing critically on g 0. Fast impurities (g < 1) be- 
have as an equivalent environment of harmonic oscilla- 
tors in dephasing the qubit, whereas for g > 1 a different 
physics emerges, dominated by memory effects, and de- 
coherence depends strongly on details of the protocols. 
We present (Fig. the decay of the qubit coherences 



T N (t) = In 



Pab(t) + Pcd (t) 



in the limit of no back-action 



Pab(0)+Pcd(0) 

of the qubit on the impurity. This is obtained by letting 
w = 0. At any fixed t = 2NAt, \T N (j)\ decreases 
monotonically when the pulse frequency 1/ At increases, 
which shows that BB effectively suppresses RTN. For 
large frequencies, 1/A£ 3> 7 (or 2A^ 3> jt), |Tjv(i)| shows 
universal behavior, scaling with g 2 . On the other hand 
for 2N 'yt qualitative differences in the behavior are 
apparent for g < 1 and g > 1. Notice that for interme- 
diate frequencies l/At < 7, the regime of experimental 
interest, BB is still able to cancel part of the noise due 
to a fast fluctuator (g < 1). For a slow fluctuator (g > 1) 
BB cancels the beats (minima in ro(£), inset of Fig.[TJ in 
the coherent dynamics 0] but besides this, it is weakly 
effective against slow RTN, despite of the semiclassical 
arguments, because there is not much to cancel. Classi- 
cal RTN causes also a systematic phase error which BB 
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FIG. 2: The purity S = ln[Tr(p Q ) 2 ] at fixed t = 87- 1 for 
protocols with N echo pair of pulses. The parameter is g — 
(fii — Qo)/j. We take e = A, v/Q = 0.2 and start from 
an eigenstate of a x . For g — 2 (slow fluctuator), S is non 
monotonic with N. Fast fluctuators (g — 0.5 and g = 0.1) 
show a more regular behavior. Inset: S scales as g 2 for AT ^> 1. 
This regime of efficient decoupling is not easily met for slow 
fluctuators (g — 2 requires At < 7 _1 /25). 



does not cancel |l9j . but can be compensated otherwise. 
Notice that the limit we discuss is the exact result for 
classical RTN, but Eq. (J5J contains also the quantum dy- 
namics of the fluctuators, including the back-action of 
the qubit. These results will be presented elsewhere [l9|. 

The physics for A ^ is even richer. We study the pu- 
rity S = lnTr(pQ) 2 , which gives deviations from unitary 
dynamics of the qubit 0. Efficient decoupling, S = 0, 
corresponds to localization in a "Zeno subspace" |5| 
which is a pure state. We study BB for generic t and 
At by diagonalizing exp(£i). Results (Fig.© show that 
for frequent pulses decoupling is achieved, S ~ 0. This 
agrees with the results of Ref. 01 for a 5 < 1 impurity. 
However decoupling slower impurities g > 1, requires 
comparatively large N. Universal behavior, S ~ g 2 is 
again found. Instead for a smaller N it may happen, 
especially for g > 1 that S is not monotonic with N, in- 
cluding the possibility that the qubit decays faster than 
in absence of pulses |3j. This is reminiscent of the inverse 
Zeno effect 13j, and it is due to the complex coupled dy- 



namics of qubit and impurity for g > 1 . 

In order to treat 1// noise we now extend our formal- 
ism to a "multi-mode" environment. We generalize the 
results at A = of Ref. |3| , to an arbitrary (non-gaussian) 
environment. The Hamiltonian is of the general form 
Eq. Q). For the evolution between two pulses at t\ and 
t 2 we use [H, <r z ] — and following the steps of Ref. |3( 
we obtain the evolution operator at t = 2N At for a BB 
protocol, S 2N (t) = [ e -<CWj.4i«r.^At e ^CW a -$,r.£}Atj* 
In the overall BB procedure a z is conserved, so we need 
only off diagonal entries of the RDM of the qubit, in the 
a z basis 

p^,(t) = p^,(Q)Tr E {s 2N (t\a)w E Sl N (t\a')} (6) 
where we assumed factorized initial conditions. Here 
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FIG. 3: BB control of 1/f noise for A = 0. The analytic 
Tjv(t) at t — 2N At (symbols - lines are guides for the eye) 
is compared with the evolution with no pulses (thick solid 
line). Noise is generated with Nfi fluctuators with rates 7; 
distributed from 10 4 Hz to 10 10 Hz. Slower fluctuators are 
ineffective 0- Noise level oc v 2 Nfi is fixed at a value typical 
of experiments in charge qubits: it is realized with coupling 
v = 9.23 ■ 10 9 Hz, for q = 6 (full triangles), 7 (circles) and 8 
(squares), with Nfi = 6 • W 17 ~ 2q scaled accordingly. Points 
for q — 6 coincide with results for gaussian noise with 1/f 
power spectrum. Crosses are stochastic Schrodinger simula- 
tions with 10 5 realizations of the 1/f process, for 5 = 7. 



S2N(t\cr) — (<j\S2n (t)\cr) expresses the conditional evo- 
lution of the environment under a well defined sequence 
of a = ±1. The trace in Eq. (jfjj factorizes if the en- 
vironment is made of noninteracting "modes" , if E are 
additive and if the initial we is factorized. If modes are 
oscillators one obtains the result of Ref. Q, which has 
been applied to a gaussian environment with 1/f spec- 
tral density . This model may have limitations |7tllO| 
in describing discrete noise sources of the solid-state, so 
we study a more realistic model, the multi-mode ver- 
sion of the Hamiltonian Eq. 0, He Jn^i an< ^ 
v d'd — > J2 V v nd\dv 13- Each "mode" is now a single im- 
purity and we take a distribution oc 7 -1 of the individual 
switching rates to produce 1/f noise 0,0- 

The contribution of each impurity to the suppression 
factor in Eq.© is calculated using the map Eq.©. The 

decay of the coherences is Fjv(t) 
Mi 



Pab(0)+ Pcd (0) 

where each plj'(t) is given by Eq.©. Results in 
Fig-EDshow that frequent pulses (curves with many sym- 
bols) suppress decoherence. Under pulsed control Tjv(t) 
changes from oc t 2 to oc t, i.e. it is described by a rate 
depending on At, as in the Zeno effect. For noise lev- 
els typical of experiments with charge qubits (Fig. © the 
pulse rate for substantial recovery is practically indepen- 
dent on v. Thus the criterion for efficient decoupling 
proposed in Ref. [2(| is not effective in this regime. The 
situation may change if a broad distribution of co upl i ngs 
is considered 0. The physics is richer for A ^ [lfllE(il | 
and compensation oil/ f noise is non monotonic for de- 
creasing At, as for a single impurity. 
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BB suppression of noise (RTN and 1//) due to quan- 
tum fluctuators is an example of general situations where 
a "structured" environment is involved. Indeed the qubit 
interacts mainly with the impurity, which is a "quantum 
filter" modulating the noise from the band. We treat this 
filter on the same footing of the qubit. 

Universal behavior in terms of the scaling parameter 
g for very frequent pulses (At -C 7 _1 ) indicates that 
when decoupling is effective details of the environment 
are unimportant. Instead in the experimentally relevant 
case of finite N (At > 7 -1 ), the different physics of slow 
(g > 1) and fast (g < 1) fluctuators manifests itself, and 
may give rise to decoupling and / or to enhancement of de- 
coherence. This picture, unexpected in the naive descrip- 
tion of BB, is reminiscent of the inverse Zeno effect 
The BB scheme we discussed prevents decoherence but 
freezes part of the dynamics. More complicated schemes 
may also allow to perform computation [2]]]. The rich 
physics we find suggests that BB may be used to ex- 
tract informations on the environment, e.g. for I// noise 
at otherwise inaccessible frequencies. Results discussed 
here are exact for classical RTN and 1// noise, but the 
formulas we presented have a broader validity: we also 
studied 0| the back-action of the qubit on the fiuctua- 
tor and 1/ / noise at general bias point, confirming the 
qualitative picture of this work. We finally stress that 
our results apply to other sources of discrete noise, as 
flux or critical current noise in flux qubits. 
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